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$E 3E. 半 干 旱地 区 矿区 的 土地 利用 格局 在 采矿 干扰 下 发 生 着 巨大 变化 ,以 全 国 八大 煤炭 生产 基地 之 一 的 山西 省 
大 同 矿 区 为 研究 对 象 ,分 析 1985 一 2015 年 土地 利用 类 型 的 时 空 变化 以 及 影响 土地 利用 变化 的 驱动 因子 ,构建 RF 


(Random Forest, RF )-FLUS (Future Land Use Simulation ,FLUS) 模 型 模拟 预测 半 干 旱 区 矿区 未 来 土地 利用 变化 ,结果 
表明 :(1) 1985 一 2015 年 ,矿区 的 林地 耕地 和 水 域 面积 减少 ,草地 和 建设 用 地 面积 增加 。(2) 林地 、 草 地 分 布 受 气候 
及 距离 水 系 和 设施 点 的 距离 影响 较 大 ;耕地 分 布 受气 候 、 高 程 及 距 水 域 . 居 民 点 的 距离 影响 较 大 ;水 域 分 布 最 重要 
的 影响 因子 是 降水 ;建设 用 地 分 布 主 要 受 生产 能 力 和 距 设 施 点 的 距离 影响 较 大 (3) FLUS 模 型 和 RF-FLUS 模 型 拟 


合 精度 均 较 高 ,但 RF-FLUS 模 型 比 FLUS 模 型 精度 更 高 ,更 接近 实际 土地 格局 变化 结果 。(4) 根据 RF-FLUS 模 型 对 


矿区 2025 年 土地 利用 变化 预测 表明 ,矿区 内 林地 .草地 条 
设 用 地 与 其 他 类 型 ( 裸 地 和 未 利用 地 ) 保 持 稳定 上 升 的 趋势 
尺度 土地 资源 优化 路 径 .促进 区 域 生态 健康 发 展 提供 有 利 的 科学 依据 。 
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地 均 呈 下 降 趋 势 ,下 降 速率 变化 不 大 ;水 域 保持 不 变 , 建 
。 本 研究 为 探究 矿区 土地 格局 复杂 动态 演变 机 制 .探索 


关键 词 : 矿区 ; 土地 利用 变化 ; 模型 预测 ;随机 森林 模型 ;驱动 因子 


日 益 加 快 的 城市 化 导致 土地 覆盖 和 土地 利用 
变化 (Land Use and Land Cover Change,LUCC ) 波 动 
明显 ” ,土地 上 覆盖 和 土地 利用 的 变化 会 造成 当地 生 
态 环境 的 破坏 ,进而 使 整个 区 域 的 生态 系统 逐渐 不 
平衡 ”。 目 前 ,对 于 土地 利用 变化 的 研究 逐步 从 现 
象 描 述 走向 机 制 解析 ,从 定性 分 析 走 向 定量 模拟 ， 
土地 利用 变化 模型 可 以 分 析 土 地 利用 变化 的 成 因 、 
探索 土地 利用 格局 演变 以 及 预测 未 来 土地 利用 的 
需求 和 分 布 5。 因 此 ,对 土地 利用 变化 动力 学 机 制 
以 及 提高 土地 利用 变化 预测 能 力 的 研究 成 为 目前 
研究 的 焦点 。 土 地 利用 变化 模型 主要 包括 基于 经 
验 统 计 原 理 的 模型 如 多 元 回归 模型 .地理 加 权 回 
归 、CLUE-S 模 型 等 “ ,基于 元 胞 自动 机 (Cellular 
Automata, CA) 原 理 的 模型 如 SLEUTH 模型 GEO- 
MOD RAISES ,基于 机 器 学 习 原 理 的 模型 如 神经 网 
络 模型 .遗传 算法 模型 等 ”, 基 于 多 智能 体 的 模型 如 
多 智能 体系 统 (Multi- Agent System, MAS) fi 789 *' 
等 。 经 验 统计 和 机 器 学 习 模 型 可 以 很 好 的 结合 于 
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感 数据 ,根据 历史 土地 利用 变化 ,建立 土地 变化 与 
时 空 要 素 之 间 的 关系 ,但 由 于 假设 驱动 力 不 变 ,所 
以 不 合适 长 时 间 段 的 预测 ,并 且 会 出 现 拟 合 过 度 问 
题 ”。 基 于 CA 原理 的 模型 采用 离散 的 空间 单元 进 
行 模拟 ,可 以 更 好 的 处 理 空间 性 信息 ,适用 于 复 
杂 地 理 空间 格局 的 模拟 和 预测 。 

以 往 CA 研究 多 是 以 逻辑 回归 为 转换 规则 进 
行 ,逻辑 回归 的 优势 在 于 运算 方便 可 以 输出 权重 便 
于 理解 ,但 精度 不 够 高 ”; 之 后 , 杨 青 生 等 ”学 者 提 
出 以 人 工 智能 和 机 器 算法 为 转换 规则 ,精度 提高 了 
但 无 法 输出 权重 。 随 着 研究 者 和 决策 者 对 土地 利 
用 分 类 未 来 格局 模拟 预测 .适宜 性 概率 获取 的 精 
度 产 生 了 更 高 的 追求 ,随机 森林 模型 结合 自身 可 以 
平衡 误差 、 复 杂 度 适中 、 分 析 驱 动因 子 权 重 等 优点 
被 广泛 应 用 于 土地 利用 研究 之 中 。 例 如 ,张大 川 
等 ”借助 随机 森林 模型 的 优势 与 CA 模型 结合 模拟 
东莞 市 土地 利用 变化 ,分 析 不 同 土地 利用 类 型 分 布 
的 各 因子 影响 重要 程度 。 马 惠 娟 等 ”在 青海 省 涅 
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水 流域 复杂 地 形 区 采用 随机 和 森林 分 类 方法 对 土地 
利用 类 型 进行 了 分 类 研究 。 陈 凯 "“ 等 利用 随机 森 
林 方 法 的 优势 结合 元 胞 自动 机 模型 通过 对 影响 城 
市 空间 分 析 的 重要 性 判别 来 研究 复杂 非 线性 的 城 
市 扩展 。 

土地 利用 变化 是 一 个 受到 自然 条 件 以 及 社会 
经 济 因 素 综合 影响 的 地 理 变化 过 程 , 由 于 矿区 的 区 
位 、 地 形 地 貌 、 开 采 方 式 和 规模 、 社 会 经 济 水 平等 
不 同 , 影 响 矿区 土地 利用 变化 的 因素 也 存在 很 大 
差异 '""。 大 同 矿区 位 于 半 干 旱 区 , 受 自 然 本 底 条 件 
的 限制 ,区 内 本 身 土地 承载 力 小 ,煤炭 资源 开发 固 
有 的 时 间 持 续 性 、 空 间 扩 展 性 和 强 干扰 性 不 仅 会 占 
用 和 破坏 大 量 土地 ,还 会 使 其 脆弱 的 生态 环境 更 
加 严重 。 如 有 段 语 凤 等 发 现 矿 区 井 工 开采 改变 了 
干旱 荒漠 区 植物 的 群落 结构 和 植被 盖 度 。 史 晓 琼 
等 研究 高 强度 采 煤 造成 的 地 下 水 位 下 降 、 地 面 沉 
降 等 生态 环境 问题 。 邻 雅 婧 等 "人 研究 发 现 黄土 高 
原 渭 北 矿区 开采 使 土壤 侵蚀 更 加 严重 。 

本 研究 基于 1985 一 2015 年 每 10 a 时 间 序 列 的 
遥感 影像 ,结合 矿区 具体 的 环境 特征 和 开采 历史 ， 
对 大 同 矿 区 近 30a 的 时 空格 局 演变 特征 及 驱动 力 
进行 分 析 , 对 未 来 2025 年 土地 利用 格局 进行 预测 。 


图 1 研究 
Fig. 1 Diagram ofthe study area 
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为 揭示 受 自 然 条 件 和 人 为 影响 的 作用 下 矿区 土地 
变化 的 规律 ,引入 随机 森林 的 方法 提高 模型 精度 更 
加 准确 预测 未 来 土地 格局 演变 ,对 于 理解 土地 利用 
格局 的 复杂 动态 演变 机 制 有 一 定理 论 意义 ,同时 对 
于 探索 矿区 土地 资源 管理 和 优化 .实施 生态 调控 策 
略 具 有 重要 实际 意义 。 


1 数据 来 源 与 方法 


1.1 研究 区 概况 

大 同 矿区 位 于 山西 省 大 同市 的 西南 方 ,地 理 位 
置 为 39°43'~40°10'N,112°31'~113°14'E, 地 形 为 向 
斜 构造 的 平缓 丘陵 ,海拔 高 度 为 1031~1964 ma W 
区 内 有 十 里 河 以 及 ”"V7" 字 形 沟谷 发 育 。 大 同 矿区 位 
于 于 旱 半 干旱 区 脆弱 地 带 , 矿 区 内 昼夜 温差 大 , 降 
水 集中 在 夏季 ,全 年 气候 干燥 ,年 均 蒸 发 量 远 大 于 
年 均 降 雨量 ,水 域 面积 小 。 矿 区 主要 含 煤 层 为 侏 罗 
系 和 石 痰 -二 和 赫 系 , 侏 罗 纪 煤层 除 西南 边 角 已 全 部 
开发 利用 。 石 痰 -二 堆 系 采 空 区 主要 分 布 在 东部 和 
南部 。 大 同 矿 区 煤 质 有 害 元 素 含量 低 ,焦油 产 率 较 
高 ,发 热量 较 高 ,是 中 国 重 要 的 优质 动力 煤 生 产 基 
地 ,开采 方式 以 井 工 开采 为 主 。 研 究 区 地 理 位置 如 
图 1 所 示 。 
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294 干旱 区 研究 3935 
1.2 数据 来 源 与 处 理 p(d, k, i) * > w,, X Sigmoid(net (d, k)) 
数据 类 型 包括 土地 利用 数据 以 及 影响 因子 数 Sax (1) 
hi -net (dk) 
据 , 影 响 因子 数据 包括 自然 因子 社会 因子 以 及 距 ? 1+e 


离 因 子 ( 表 1)。 土 地 利用 数据 是 由 30 m 分 辩 率 的 
Landsat 系列 遥感 影像 数据 解 译 获 取 , 并 按照 国家 标 
准 重 分 类 为 六 类 。 高 程 数 据 由 中 国 科 学 院 资源 环 
境 科 学 数据 中 心 (http:/www.resdc.cn) 获 取 ,并 从 中 
提取 出 坡 向 和 坡度 。 气 温和 降水 数据 由 中 国 气象 
数据 网 (http://data.cma.cn) 气 温和 降水 站 点 数据 进 
行 空 间 持 值 获得 。 煤 矿 生 产能 力 由 实地 调查 获 
得 。 人 口 数据 在 Ye 等 研究 基础 上 将 精度 提升 到 
30 m。 距 离 因子 通过 采用 ArcGIS 10.5 中 欧 氏 距离 
算法 计算 得 出 ,POI 数 据 采用 网 络 爬 虫 技 术 获 得 。 
1.3 研究 方法 
13.4 FLUS 模型 FLUS 模 型 耦合 了 系统 动力 学 模 
型 和 元 胞 自动 机 模型 ,并 在 元 胞 自动 机 模型 原 有 基 
础 上 引入 人 工 神 经 网 络 (Artificial Neural Network , 
ANN) 构 建 转换 规则 ,具有 耦合 人 类 活动 与 自然 效 
应 的 各 类 土地 类 型 之 间 相 互 苋 争 和 转化 的 模拟 优 
势 ””。FLUS 模 型 包括 如 下 两 个 部 分 : 

第 一 部 分 是 通过 ANN 模块 训练 样本 数据 并 构 
建 本 区 域 的 土地 利用 类 型 空间 发 生 概 率 。 人 工 神 
经 网 络 的 优势 在 于 能 够 通过 多 次 迭代 学 习 和 拟 合 
输入 数据 与 训练 目标 之 间 的 复杂 关系 中。 在 ANN 
的 输入 层 中 ,神经 元 与 输入 变量 相对 应 ;输出 层 中 
的 神经 元 代表 特定 地 类 在 网 格 单元 内 的 发 生 概率 ， 
在 特定 单元 格 内 该 值 越 高 表示 目标 地 类 的 发 生 概 
率 越 高 ”。 根 据 公 式 可 得 出 ,在 时 刻 第 i 种 地 类 在 
WHE d EEE WER: 


式 中 : p(d. k, i) RIRES k Z EIEEE d E i TARS 
适应 性 概率 ; w, 表示 输出 层 与 隐藏 层 间 的 自 适 应 
权重 ; nei,(d, i) 表示 隐藏 层 中 神经 元 j 接 收 到 的 信 
号 ; Sigmoid 为 输入 层 与 隐藏 层 之 间 的 连接 激活 

第 二 部 分 是 基于 轮 盘 赌 选 择机 制 的 自 适应 惯 
性 竞争 机 制 融 合 CA 的 空间 运算 功能 ,调整 宏观 土 
地 需求 与 当前 现 有 土地 之 间 的 数量 并 达到 目标 
值 , 解 决 了 以 往 模 型 在 处 理 复杂 地 类 间 的 相互 影 
Wo. XA: 


Intertia; ', 


ifs |< |s 


Intertiaj ' X i if|5; | «|S;"| «0 (2) 


Intertia; = 


Intertia; ! X 


S : i-2 i-l 
ge if0«|S;^|«|S;"| 


式 中 : Intertia; 为 i 时 刻 第 种 地 类 的 自 适应 惯性 系 
Jo Sc Si 代表 广 1 -2 时 刻 第 天 种 地 类 的 需求 与 
当前 土地 数量 的 差 值 。 转 移 成 本 由 土地 利用 历史 
数据 与 经 验 得 出 ,表示 为 se, sse 
综 上 可 以 将 综合 适宜 性 概率 计算 公式 表示 为 : 

TP;,, - p(d, i, k) X OQ x Intertia, X (1 —sc,, ,,) (3) 
1.3.0 RF-FLUS 模 型 ”随机 森林 模型 是 一 个 在 数据 
挖掘 能 力 ,转换 规则 准确 性 及 其 稳定 性 方面 具有 自 
吴 优 势 的 利用 多 个 决策 树 生 成 的 集成 智能 算法 。 
只 需要 对 给 定 样 本 进行 训练 ,克服 神经 网 络 过 度 拟 
合 的 问题 人 *。 采 用 RF 代 蔡 ANN 在 FLUS 模 型 内 的 转 
换 规则 ,构造 24 棵 决策 树 ,变量 数 为 4, 每 个 样本 抽取 
2/3 的 数据 作为 袋 内 数据 用 来 训练 模型 ,剩余 1/3 的 


表 1 数据 来 源 


Tab.1 Data source 


数据 类 型 数据 名 称 数据 来 源 
自然 因子 高 程 中 国 科学 院 资源 环境 科学 数据 中 心 (http://www.resdc.cn ) 
坡度 、 坡 向 基于 高 程 (DEM ) 数 据 提取 
气温 、 降 水 中 国 气象 数据 网 (http://data.cma.cn) 
社会 经 济 因 子 煤矿 生产 能 力 采用 实际 调查 获取 
AH 中 国 科学 院 资源 环境 科学 数据 中 心 (http://www.resde.cn) 
距离 因子 距 公 路 \ 铁 路 .水系 距 离 KH ArcGIS 10.5 中 欧 氏 距离 算法 计算 
距 设 施 点 、 居 民 点 距离 POI 数 据 采 用 网 络 怜 虫 技术 获取 服务 设施 
土地 利用 数据 1985 年 .1995 年 .2005 年 .2015 年 遥感 影像 中 国 科学 院 资源 环境 科学 数据 中 心 (http:/www.resdc.cn ) 
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数据 用 作 袋 外 数据 来 估算 误差 并 验证 决策 树 的 精 
度 , 袋 外 数据 越 小 精度 越 高 。 样 本 和 变量 的 选择 可 
以 降低 各 决策 树 之 间 相 关 性 ,为 之 后 更 准确 的 对 大 
同 矿区 土地 利用 变化 进行 模拟 和 预测 提供 基础 。 

1.3.3 平均 基尼 系数 下 降 法 ”选用 平均 基尼 系数 下 
降 法 来 计算 并 分 析 影 响 大 同 矿区 土地 利用 格局 不 
同 因素 的 影响 程度 。 平 均 基 尼 系 数 下 降 法 是 对 每 
颗 决 策 树 按照 基尼 的 不 纯度 (Gini Impurity ) 对 特征 
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1995 一 2005 年 矿区 走向 成 熟 阶段 ,生产 过 程 逐 
渐 机 械 化 ,开采 规模 进一步 扩大 。 耕 地 主要 零散 分 
布 于 矿区 边缘 ,该 时 期 耕地 面积 减少 约 20% ; 随 着 
退耕 还 林 还 草 政策 的 推行 ,草地 成 为 矿区 所 占 比 例 
最 大 的 土地 利用 类 型 , 占 总 面积 的 60% ;水 域 面 积 
急剧 下 降 , 由 1.04% 减 少 至 0.1% ;建设 用 地 减少 
0.67% ,其 他 类 型 增加 0.35%。 

2005 一 2015 年 矿区 处 于 转型 期 ,政府 政策 兵 持 
生产 全 面 改 造 , 矿 井 生产 逐渐 规模 化 ,同时 矿区 城 


进行 排序 ,特征 的 重要 性 直接 反映 了 该 特征 参与 模 
型 计算 过 程 中 的 重要 性 。 每 个 特征 的 重要 性 为 0~1 
之 间 的 浮 点 数 ,数值 越 大 表明 该 特征 对 于 整个 模型 
越 重要 。 通 过 遍历 全 部 树 节 点 ,计算 选 入 的 每 个 特 
征 变量 的 基尼 系数 下 降 总 和 表示 该 变量 对 最 终 格 
局 的 影响 程度 ? 。 

1.3.4 模型 精度 验证 ”模型 精度 验证 采用 混 清和 矩阵 
Kappa 系 数 检验 所 定量 分 析 验 证 模型 精度 。 混 消 拢 
阵 对 角 线 百分比 越 高 ,代表 精度 越 高 。 


2 结果 与 分 析 


2.1 土地 利用 格局 时 空 变 化 特征 

由 1985 一 2015 年 大 同 矿 区 土地 利用 格局 空间 
变化 和 面积 变化 图 可 知 (图 2, 图 3),1985 一 1995 年 
矿区 处 于 成 长 阶段 ,耕地 面积 增加 了 40.75 km; 草 
地 广泛 分 布 于 研究 区 内 ,变化 率 为 -2.87% ;林地 主 
要 分 布 在 中 部 以 及 东部 边缘 ,面积 增加 3.37 km? , 林 
草地 数量 仅 次 于 耕地 面积 ;水 域 主要 分 布 在 研究 区 
西北 部 的 十 里 河 , 呈 带 状 分 布 水 量 较 多 ,面积 增加 
了 3.53 km2; 建 设 用 地 和 其 他 类 型 多 分 布 在 北部 地 
势 平坦 处 ,面积 增加 了 0.27% 。 


(a) 1985 年 (b) 1995 年 


(c) 2005 年 


市 工业 化 程度 不 断 加 快 , 因 此 该 期 土地 利用 特点 主 
要 表现 为 建设 用 地 面积 大 幅 增 加 ,由 3.09% 增 加 到 
7.09%。 耕 地 、 林 地 草地 以 及 水 域 地 类 由 于 人 类 活 
动 的 影响 面积 都 有 所 下 降 ,其 他 类 型 面积 略微 上 
升 ,增加 0.29%。 
2.2 不 同 历史 时 期 土地 利用 转移 分 析 

H 1985—2015 年 大 同 矿区 的 土地 利用 转移 可 
知 ( 表 2) , 1985—1995 年 由 于 采矿 活动 和 人 类 生活 
的 影响 ,土地 利用 转移 主要 为 水 域 大 面积 转变 为 林 
地 、 草 地 和 耕地 ,建设 用 地 有 6.72% 转 向 林地 ; 
1995—2005 年 矿区 开始 大 规模 开采 , 毁 林 采矿 JE 
筑 公 路 实现 晋 煤 外 运 ,大 量 占用 耕地 和 林地 ,土地 
利用 的 变化 主要 表现 为 各 土地 利用 类 型 向 草地 的 
转移 ,耕地 向 草地 转移 高 达 72.58%, 林 地 向 草地 发 
生 转 移 55.93% ,水 域 建设 用 地 和 其 他 用 地 向 草地 
转移 均 为 30% 以 上 。 水 域 面 积 向 耕地 转移 较 多 , 达 
33.39% ; 2005—2015 年 土地 利用 类 型 转移 主要 表现 
为 水 域 转 为 其 他 类 型 ,转移 94.35% ,少量 耕地 和 水 
域 转 为 曹 地。 其余 土 地 利用 类 型 几乎 没有 发 生 转 
移 。30 a 间 大 同 矿区 土地 利用 的 变化 趋势 是 各 土地 
利用 类 型 癌 草 地 转移 ,其 中 ,林地 向 草地 转移 最 多 ， 


(d) 2015 年 


图 例 
mm PAM mm EUR 四 耕地 


N 
A 0 10km mki 四 其 他 中 建设 用 地 


图 2 1985 一 2015 年 大 同 矿区 土地 利用 空间 分 布 


Fig.2 Spatial distribution of land use in Datong coal mining area from 1985 to 2015 
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图 3 大 同 矿区 不 同时 期 土地 类 型 面积 变化 
Fig. 3 Changes of land type area in Datong coal mining 


area in each year 


其 次 是 耕地 。 
2.3 土地 利用 变化 的 影响 因素 分 析 

根据 土地 利用 变量 的 重要 性 ,进一步 分 析 大 同 
煤矿 区 土地 利用 变化 的 影响 因子 (图 4) , 林 草 地 的 
分 布 主要 受到 气候 ` 距 水 系 和 设施 点 距离 影响 ; 耕 
地 分 布 受到 矿区 DEM 气候 、 距 水 系 \、 设 施 点 以 及 居 
民 点 的 距离 等 因素 影响 。 这 主要 由 于 干旱 区 气候 


耕种 将 耕地 分 布 在 距离 居民 点 较 近 的 地 方 ;水 域 最 
重要 的 影响 因子 是 降水 因子 ,降水 条 件 直接 决定 水 
域 分 布 范 围 ; 建 设 用 地 很 大 程度 上 受到 生产 能 力 和 
距 设 施 点 距离 的 影响 ,矿区 煤炭 生产 能 力 越 高 ,各 
种 厂房 . 排 土 道路 .煤矿 职工 办 公 住 宿 区 均 相 应 增 
加 。 餐 馆 、 学 校医 院 、 休 闲 娱乐 场所 等 设施 点 的 距 
离 也 直接 决定 建设 用 地 分 布 。 其 他 类 型 分 布 受到 
距 水 系 距离 的 影响 ,水 系 周围 的 用 地 极 易 被 改变 为 
其 他 类 型 。 
2.4 土地 利用 格局 模拟 及 未 来 预测 

分 别 采用 FLUS 和 RF-FLUS 模 型 ,对 大 同 矿区 
土地 利用 变化 进行 建 模 , 采 用 Kappa 系数 对 比 精度 
结果 。 首 先 对 2015 年 大 同 煤 矿区 土地 利用 格局 进 
行 拟 合 ,结果 表明 ,FLUS 模 型 和 RF-FLUS 模 型 模拟 
得 到 的 2015 年 矿区 土地 利用 格局 与 2015 年 影像 解 
译 得 到 的 实际 土地 利用 格局 相 比 (图 5a~ 图 5c), 拟 
合 程度 均 较 高 。 将 模拟 结果 与 实际 土地 利用 做 混 
消 和 矩阵 精度 验证 ( 表 3)。FLUS 模型 Kappa 系 数 为 
0.92, OA 系数 为 0.95; RF-FLUS 模型 Kappa 系数 为 
0.93 ,OA 系数 为 0.96 ,就 模拟 精度 而 言 ,RF-FLUS Bé 
型 拟 合 效 果 更 好 ,更 加 符合 实际 土地 利用 格局 。 

以 2015 年 为 基准 年 ,采用 RF-FLUS 模 型 模拟 矿 


干燥 ,耕地 需要 便利 的 灌溉 条 件 ,而 且 农民 为 方便 


区 2025 年 的 土地 利用 空间 分 布 结果 可 知 ( 图 6)。 


表 2 1985—2015 年 大 同 矿 区 土地 利用 转移 概率 矩阵 


Tab.2 Land use transfer probability matrix in Datong coal mining area from 1985 to 2015 /% 
年 份 土地 利用 类 型 林地 草地 耕地 水 域 建设 用 地 其 他 
1985 一 1995 林地 86.88 3.69 8.54 0.45 0.44 一 
草地 5.26 77.63 16.49 0.31 0.20 0.11 
耕地 14.87 8.82 73.60 0.58 1.94 0.19 
水 域 17.31 23.64 5.41 52.47 1.17 一 
建设 用 地 6.72 3.18 6.27 4.02 77.55 2.25 
其 他 26.68 1.64 一 一 2.35 69.33 
1995—2005 林地 10.44 55.93 312 0.10 2.07 0.27 
草地 35.7 49.14 13.26 0.14 1.07 0.68 
耕地 13.36 72.58 11.59 0.12 1.75 0.60 
水 域 21.23 3128 33.39 0.33 13.47 0.30 
建设 用 地 4.17 32.89 27.64 0.05 32.92 2.33 
其 他 7.02 47.12 7.08 0.37 28.75 9.66 
2005 一 2015 林地 99.20 一 0.06 一 0.30 0.44 
草地 = 97.54 0.05 一 2.24 0.18 
耕地 - 1.03 85.99 - 12.97 0.01 
水 域 - 2.12 3:52 - 94.35 
建设 用 地 一 一 - - 100.00 - 
其 他 = E - - - 100.00 
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图 4 2015 年 大 同 矿 区 土地 利用 类 型 影响 因子 的 重要 性 


Fig.4 The importance of different types of influencing factors in Datong coal mining area in 2015 


(a) 实际 土地 利用 格局 I (b) FLUS 模 型 模拟 结果 7 (c) RF-FLUS 模 型 模拟 结果 7 


0O lOkm 图 例 meg 林地 ENSE 草地 [9 耕地 ok ll 建设 用 地 E 其 他 


图 5 2015 年 大 同 矿区 土地 利用 格局 与 模型 模拟 结 


Fig. 5 Land use pattern and model simulation results in Datong coal mining area in 2015 
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表 3 RF-FLUS 模型 与 FLUS 模型 混淆 矩阵 
Tab.3 Confusion matrix of RF-FLUS model and FLUS model /% 
RF-FLUS 模 型 FLUS 模 型 
林地 “草地 He 水 域 建设 用 地 ”其 他 林地 草地 “耕地 水域“ 建设 用 地 ”其 他 
林地 95.03 1.73 1.35 0.00 0.00 0.20 ”林地 93.25 2.14 1.80 0.00 0.00 0.33 
草地 2.35 9571 2.60 0.00 0.00 571 ”草地 2.904 . 95.33 2.32 10.00 0.00 6.03 
耕地 2.02 120 94.57 0.00 0.00 026 ”耕地 2.54 1.62 93.52 0.00 0.00 0.26 
水 域 0.00 0.00 0.00 — 100.00 0.00 0.00 水域 0.00 0.00 0.00 90.00 0.00 0.00 
建设 用 地 。 028 1.17 1.40 0.00 100.00 013 ”建设 用 地 0.66 0.85 2.05 0.00 100.00 0.20 
其 他 0.32 0.19 0.08 0.00 0.00 9370 其 他 0.61 0.06 0.31 0.00 0.00 93.18 
总 计 100.00 100.00 100.00 100.00 100.00 100.00 总 计 100.00 100.00 100.00 100.00 100.00 100.00 
(a) 2015 年 (b) 2025 年 的 同时 造成 对 土地 的 破坏 ,从 而 改变 了 土地 利用 


图 例 
N 林地 m Eib 耕地 
A 0 13km DO m RUD m C 


图 6 大 同 矿区 2015 年 和 2025 年 土地 利用 分 布 格局 
Fig. 6 Distribution pattern of land use in Datong coal mining 
area in 2015 and 2025 


2025 年 大 同 煤矿 区 林地 .草地 和 耕地 均 呈 下 降 态 
势 , 林 地、 草地 和 耕地 面积 分 别 下 降 了 3.01 km, 
14.97 km 和 12.7 km? ,水 域 面积 基本 保持 不 变 ,建设 
用 地 与 其 他 类 型 面积 稳 态 上 升 ,建设 用 地 由 2015 年 
的 130.75 km? 增加 到 2025 年 的 157 km ,增加 了 
26.26 km? ,其 他 类 型 面积 预计 增加 44.12 km 
3 结论 

本 文通 过 研究 大 同 矿 区 不 同 阶段 土地 利用 格 
局 的 时 空 变化 ,进一步 探讨 不 同 土地 利用 模型 的 模 
拟 精 度 优势 以 及 影响 土地 利用 格局 变化 的 主要 驱 
动因 素 ,并 进行 了 预测 分 析 , 得 到 如 下 结论 : 

(1) 大 同 矿 区 侏 罗 系 煤层 已 基本 全 部 开发 利 
用 ,石炭 -二 又 系 煤层 开采 范围 主要 在 北部 和 南部 ， 
受 采 矿 影 响 ,1985 一 2015 年 大 同 矿 区 土地 利用 格局 
发 生 剧 烈 转 变 , 主 要 表现 为 草地 建设 用 地 以 及 其 
他 类 型 增加 ,林地 和 耕地 先 增加 后 减少 ,以 减少 为 
主 ,水 域 面积 基本 被 侵占 。 矿 区 土地 利用 变化 是 以 
资源 开采 为 原动力 的 动态 时 空 演变 过 程 , 采 掘 资源 


结构 。 

(2) 矿区 内 不 同 土地 利用 类 型 的 驱动 因素 不 
同 。 林 地 、 草 地 主要 受气 候 有 影响 以 及 距 水 系 和 设施 
点 的 距离 影响 ;耕地 分 布 受到 高 程 .气候 及 距 水 系 、 
设施 点 .居民 点 的 距离 等 多 因素 综合 影响 ;水 域 最 
重要 的 驱动 因素 是 降水 。 建 设 用 地 受到 生产 能 
和 上 距 设施 点 的 距离 的 影响 。 其 他 类 型 主要 受 距 水 
系 距 离 的 影响 。 矿 产 资 源 开发 是 土地 利用 变化 的 
主要 驱动 力 , 此 外 气候 地 形 等 自然 因素 也 影响 着 土 
地 利用 变化 。 

(3) RF-FLUS 模 型 的 应 用 ,可 以 更 加 准确 的 预 
测 未 来 土地 利用 格局 。 相 比 于 FLUS 模 型 中 ANN 模 
块 出 现 过度 拟 合 的 现象 ,RF-FLUS 模型 能 发 挥 随机 
森林 算法 精度 高 的 优势 ,提高 预测 未 来 土地 利用 格 
局 变化 的 能 力 , 对 预测 小 尺度 区 域 的 土地 利用 格局 
具有 重要 借鉴 意义 。 
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Abstract: Land use and land cover change (LUCC) is a critical factor contributing to regional land planning and 
ecological environmental protection. LUCCSs associated with human activities, such as coal mining, have 
increased the tension in human-environment interactions. This study explored the spatio-temporal LUCC and its 
driving factors from 1985 to 2015 in the Datong mining area, which is one of the eight major coal production 
bases in China. Furthermore, the random forests (RF)-future land use simulation (FLUS) model was proposed to 
explicitly simulate the spatial trajectories of LUCCSs for the year 2025. The results showed that (1) forestland, 
cropland, and watershed areas kept decreasing, while grasslands and construction lands kept increasing from 
1985 to 2015. (2) Climate, elevation, and distance were the most influential factors for the distribution of 
croplands, forestlands, and grasslands, while precipitation was the most important factor for the distribution of 
watersheds. Coal production capacity and distance from facilities were the most influential factors for the 
distribution of construction lands. (3) Both the FLUS and RF-FLUS models had a good fitting accuracy, whereas 
the RF-FLUS models had a satisfactory kappa index and OA index. (4) Land use simulations for the year 2025 
based on the RF-FLUS model indicated that croplands, forestlands, and grasslands will decrease, while 
construction lands will increase. This study provides a scientific reference for understanding the complex and 
dynamic evolution of LUCC, exploring the possibilities for the way of land resources, and promoting the 
sustainable development of coal mining areas. 
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